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Abstract 

We derive relativistic hydrodynamics from quantum field theories by assuming that the den¬ 
sity operator is given by a local Gibbs distribution at initial time. We decompose the energy- 
momentum tensor and particle current into nondissipative and dissipative parts, and analyze their 
time evolution in detail. Performing the path-integral formulation of the local Gibbs distribution, 
we microscopically derive the generating functional for the nondissipative hydrodynamics. We also 
construct a basis to study dissipative corrections. In particular, we derive the first-order dissipative 
hydrodynamic equations without a choice of frame such as the Landau-Lifshitz or Eckart frame. 

PACS numbers: 47.75.-|-f, 47.10.-g 
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I. INTRODUCTION AND SUMMARY 


Hydrodynamics universally describes the spacetime evolution of charge densities of sys¬ 
tems such as energy, momentum, and particle number [1]. It does not depend on microscopic 
details of systems, whose application covers branches of physics from condensed matter to 
high-energy physics. Among them is illuminating the recent success of relativistic hydrody¬ 
namics in describing the evolution of the quark-gluon plasma (QGP) created in heavy-ion 
collision experiments [2-8]. 

The hrst-order relativistic hydrodynamic equations, that is, the relativistic version of the 
Navier-Stokes equations, which suffer from the violation of causality, have been derived by 
Eckart [9] and by Landau and Lifshitz [1]. The second-order equations, which resolve the 
causality problem by introducing a hnite relaxation time, were derived hrst by Muller [10] 
and also by Israel and Stewart [11]. After the aforementioned success of relativistic hydro¬ 
dynamics in describing the QGP, a lot of work concerning the derivation of hydrodynamic 
equations has been progressively carried out, in which the hydrodynamic equations are for¬ 
mulated based on the kinetic theory [12-19], the fluid/gravity correspondence [20-23], the 
phenomenological extension of the nonequilibrium thermodynamics [24, 25], and the pro¬ 
jection operator method [26, 27]. Also, a signihcant method has recently been developed 
in which the equilibrium-generating functional for the nondissipative hydrodynamics is con¬ 
structed only by respecting symmetries of systems [28, 29]. 

The aim of this work is to derive the dissipative relativistic hydrodynamic equations 
from quantum held theories. Our approach is based on the recent development of the 
nonequilibrium statistical mechanics [30], which is essentially equivalent to the nonequi- 
librium statistical operator method [31, 32]. By performing the path-integral formulation 
of the Massieu-Planck functional, we present the hrst microscopic justihcation of the gen¬ 
erating functional method [28, 29] for nondissipative parts. This enables us to justify a 
generalized argument by Luttinger [33] , in which the spatial distribution of the temperature 
is interpreted as an auxiliary external gravitational potential. We also study the dissipative 
corrections to relativistic hydrodynamic equations by using our method. Although we re¬ 
strict ourselves to hrst-order equations in this paper, our formulation provides a solid basis 
to proceed to the higher orders in the derivative expansion. 

In the rest of this section, we briehy summarize our result. The relativistic hydrodynamic 
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equations are based on the continuity equations: 


( 1 ) 

( 2 ) 


= 0 , 

= 0 . 

Here is the covariant derivative. T'^'' and are the energy-momentum tensor and 
particle current, respectively. They are decomposed into nondissipative and dissipative 
parts 


_j_ 


r = 4 + 


( 3 ) 

( 4 ) 


In the leading order of the derivative expansion, the nondissipative terms have the form of 
a perfect fluid: Tq'" = (e + p)u^u’^ + and Jq = nu^. Here e denotes the energy density, 
p the pressure, n the particle density, and the fluid four-velocity. and 5J^ represent 

the dissipative parts. In our formalism, the dissipative terms are given as 

( 5 ) 

5J^ = ( 6 ) 

in the leading order of the derivative expansion , where (3^ = l3u^ with the inverse tem¬ 
perature /3, and u = (3p with the chemical potential p. Here, (, p, and n denote the bulk 
viscosity, the shear viscosity, and the diffusion constant, respectively, whose microscopic 
expressions are given by the Kubo formulas, Eqs. (83) — (85). We introduced the spatial 
projection operator Pj^ = 5^ + v^riy, and , where denotes the normal vec¬ 

tor for an isochronous hypersurface and the time vector with = — 1. These spatial 
projection operators satisfy P^^v^ = P^^n^ = 0 and PpP^ = . We also dehned tensors 

with angle brackets as the traceless symmetric projected parts, which are given explicitly as 

= Ip;pS(^<, 9 + As„) - (7) 

where d is the spacetime dimension, is a symmetric tensor and satishes = Pj^. 

We emphasize here that the above constitutive relations, Eqs. (5) and (6), with the Kubo 
formulas, Eqs. (83) — (85), are obtained without choosing any frame; this is an advantage 
of our new formulation. The particular choice of and reproduces the dissipative 
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hydrodynamic equations in the known frame. For example, we reproduce the Landau- 
Lifshitz frame if we choose . 

This paper is organized as follows: In Sec. II, we review the local thermodynamics. 
In Sec. IIC, we derive the path-integral formulation of the Massieu-Planck functional on a 
hypersurface. In Sec. Ill, we discuss the time evolution of hydrodynamic variables, and derive 
self-consistent equations giving constitutive relations. In Sec. IV, we discuss the derivative 
expansion of the hydrodynamic equations in a frame-independent way. Section. V is devoted 
to a discussion. 


II. LOCAL THERMODYNAMICS ON A HYPERSURFACE 

In this section we discuss the local thermodynamics on a spacelike hypersurface in order 
to construct relativistic hydrodynamic equations in a covariant way. In Sec. II A, we hrst 
summarize geometric aspects of the spatial hypersurface used in this paper. In Sec. IIB, 
we introduce several concepts such as the local Gibbs distribution and the entropy current 
operator based on Refs. [31, 32, 34, 35]. In Sec. IIC, we derive the path-integral formulation 
of the Massieu-Planck functional on the hypersurface. The Lagrangian is written as that 
in the curved spacetime background helds, whose metric consists of the local temperature 
and the fluid four-velocity. We show that the metric has Kaluza-Klein gauge symmetry in 
addition to {d — l)-dimensional diffeomorphism invariance [28]. 



FIG. 1. Illustration of the Arnowitt-Deser-Misner (ADM) decomposition of the spacetime. Sj 
denotes a spacelike hypersurface parametrized by i{x) = const, is a vector normal to the 
hypersurface. Introducing the lapse function N{x) and the shift vector N^{x), we decompose the 
time vector as d^x^ = Nv^ = Nn^ + N^. 


4 






A. Geometric preliminaries 


As a technical preparation, we first summarize the geometric aspects of spacelike hyper¬ 
surface in this subsection. Let us consider spatial slicings on a general curved spacetime 
with a metric gfj,„ and parametrize the spacelike hypersurface by i. We also introduce the 
spatial coordinates x on the hypersurface. In other words, we dehne a spacelike hypersurface 
Sf by the i{x) = const, surface, and introduce spatial coordinates x = x{x), where x is a 
general coordinate (see Fig. 1). To discuss dynamics on such a spacelike hypersurface, it is 
convenient to introduce a timelike unit vector as 

n^{x) = —N{x)d^t{x) with N{x) = {—dH{x)d^t{x))~^^‘^. (8) 


Here we normalize as = — 1 and is future oriented. A > 0 is the lapse function. 
We use the mostly plus convention of the metric, e.g., the Minkowski metric is = 
diag(—1,1,1, ■ ■ • , 1). The induced metric on the spacelike hypersurface is then 

Ifiu = QiJiv + . (9) 

We also introduce the shift vector by the decomposition 


diX^{t, x) = Nn^ + with = 0 . 

In the coordinate system {i,x), and are given explicitly by 

np = {-N, 0 ) , yoi = Tio = 9oi = 9lo. lij = 9ij, = 


( 10 ) 


0 

N^gOi 


(11) 


The metric gpp takes the form of the Arnowitt-Deser-Misner (ADM) metric. 


dx^ dx'^ I —N‘^ + WA* W 


9p.p 9^,^ 


dxi^ dx'^ 


Ni 


9^’^ = 


lij 


_iV-2 iV-2W 


( 12 ) 


Here = TfjNK y*l is the inverse of jij and satishes = S^. The d-dimensional 

volume element is given by 


J d^x\/^ = j d^xN ^ with y = dety^j , 
whereas the volume element on the spacelike hypersurface Sf is 

f d'Li= [ d'^x\/^^5(t — t{x))N~^{x) = f d'^~^Xy/^. 


(13) 


(14) 
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It is also convenient to introduce a vector proportional to Eq. (10) as 

x) with = —1. (15) 

Using and we define a spatial projection operator as 

Pj; = 6^Pv>^n^ with = 0, = 0, P^PP = Pj^. (16) 

Its concrete form in the coordinate system (t, x) is given by Pp = diag(0,1,1,..., 1). We 
will use this projection operator in Sec. III. We note that such an operator often appears in 
the context of Newton-Cartan geometry (see, e.g., Refs. [36, 37] ^). 

B. Local Gibbs distribution 

We next introduce a density operator representing a local thermal equilibrium state, 
and review the thermodynamics on the hypersurface [31, 32, 34, 35]. We start with global 
thermal equilibrium on the Minkowski space, in which the density operator for an arbitrary 
inertial frame of reference is given as the Gibbs distribution, 

where parameters are (3^ = with the inverse temperature f3, the fluid four-velocity of the 
system normalized by = — 1, and u = (3p with the chemical potential p. P^ and N 
denote energy-momentum and number operators, respectively. The Massieu-Planck function 
4/(/9^,i/) = Intr exp[/3^P^ -|- uN] determines the normalization of the density operator peq. 
At the rest frame of medium, = (1, 0 ), and thus peq(/d^, u) = exp [—/3(P — pN) — 4/(/3, z/)] 
are satisfied. 

We then generalize the global Gibbs distribution (17) to a local form in a coordinate- 
invariant way. For this purpose, let us consider thermodynamics on the spacelike hypersur¬ 
face, Sp introduced in the previous subsection. For generality, we leave the metric of 
the spacetime as a general curved one. On the hypersurface, we introduce a local Gibbs 
distribution plg[^;A] as 

Plg[UA] = exp(-5[f; A]) with ^[t; A] = P[f; A]T[t; A], (18) 

^ Our normalization = — 1 has the opposite sign compared to that in Refs. [36, 37] 
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where K[t]\] is defined by 

k[i-\] = - I dE,^\kx)jnk = - I dE,, + u{x)r{x)). (19) 

Here we introduced dT,t^ = —dTiin^. A“ and denote sets of parameters, A“(a;) = 
{/9^(x), and of current operators, Ja{.x) = {T'^(a;), J^(a;)}, respectively. Just as in 

the global case (17), the Massieu-Planck functional J/[t; A] = lntrexp(—iP[t; A]) determines 
the normalization of the density operator pi^Q. For constant parameters and = (—1,0), 
the local Gibbs distribution reproduces the global one (17). We note that the definition here 
is coordinate invariant by construction. 

The charge density operators on the hypersurface, Ca(x) = {p^(x), ?7'(x)}, read p^x) = 
—ny{x)T''kx) andh'(x) = —nu{x)J^{x). Their expectation values, (ca(x))j'^ = tr [/3 lg[^; '^]ca(a^)], 
are obtained from the variation of \h[7; A] with respect to A“(x) on Sf, 

Caix) = (Ca(x))p = ^]- (20) 

The entropy is defined by 


S[t;c\ = -trpLG[t; A] lnpLG[^;A] 

= {S[t;X])^^ ( 21 ) 

= -J dSfA“c, + T[7;A]. 


The entropy is a functional of c^, not A“, which can be confirmed by conducting the variation 
of S with the fixed i, 


6S 




JT[t;A] 

JA“ 



j dEiX^Sc,. 


( 22 ) 


The parameters are obtained as 


(23) 


For later purposes, we introduce -0^ such that 


T[7;A] 


/ 


= 


dEiij, 


(24) 
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where 'ifj = which satishes 


dip = CadX°' = p^dfd^ + n'dv, (25) 

up to the covariant total derivative that does not contribute to 5\1/. As will be seen in 
Sec. IV A 1, in the leading order of derivative expansion, we can write ip^ as ip^ = P'^p{f3, v) 
with the pressure p. We note that there is an ambiguity in the dehnition of ip^ because d' 
is invariant under the transformation -0^ —>■ -0^ + with a function satisfying n^g^ = 0. 
Introducing the entropy current operator, 

s^^ = -X^jpp + lP^^ = -/3'^f^p - i/> + (26) 

the entropy reads 

S = j dT.i^s^ = j dEtS, (27) 

where and s = —n^s^ = —X°'Ca + ip = —fd^Pfi — I'n' + ip. The entropy density 

s satishes the thermodynamic relation, ds = —X^'dca = —— udn', up to the covariant 
total derivative. The divergence of the entropy current is 


( 28 ) 

where we used the continuity equations = 0. In order to evaluate V^ip^, let us 

consider the derivative of \h[f. A] with respect to t, which reads 


= {di I 

= (y rfEtivvqA“j:))7 

= j dE,-iV(V,.A“)(j:)p, 

where we again used the continuity equations. We also used 


(29) 


Sf / = 


dT.iNVJ‘‘, 


( 30 ) 


for an arbitrary smooth function f^{x) (see Appendix A). 


From Eq. (29), we obtain the divergence of "0^ as 

v„r = (v^A-x (31) 

Then, the divergence of the entropy cnrrent operator reads 

- (V;.z/)<5>, (32) 

where 50 = O — {O)^^. The entropy prodnction rate (V^s^) is in general nonzero. When 
we decompose the expectation value of the current as {Ja) = 

can be identihed as the nondissipative part because it does not contribute to the entropy 
production rate, while {557^) can be identihed as the dissipative part. 


C. Path-integral formulation of Massieu-Planck functional and thermal metric 


In this subsection, we derive the path-integral formula for the Massieu-Planck functional 
T. We show that the action has a form in the curved spacetime background, whose metric 
depends on parameters and u. We also show that the result is in accordance with those of 
recent studies, in which the Massieu-Planck functional is derived on the basis of symmetric 
and scaling properties [28, 29]. Although we only consider a neutral scalar held here, the 
discussion covers the essential feature of the Massieu-Planck functional. 

In the coordinate system {t, x) with the ADM metric (12), the Lagrangian for a neutral 
scalar held 0 reads 
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2JV2 - \di4,d-,4: - V{4>), 


(33) 


where V{(j)) denotes the potential term. The canonical momentum 7i{x) is 7t{x) = 
—g^^dp(j){x) = N~‘^{dQ(j) — N^dicp), which satishes the canonical commutation relation, 
[0(S), 7r(S')] = i6{x — x')/{Ny/j). We obtain the energy-momentum tensors as 


N' 
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T'l = T^dicp - c = —TT^ + N'-ndicj) + —di^id-jcl) + D(0), 

t'i = 

By using the standard technique of the path integral, we have 
tre“'^= [ d(j){(j)\e~^\(j)) 


(34) 

(35) 


= / P0T>7r exp 


'•00 


dT[i I d'^ ^xN^dT-(j){T,x)'K{T,x) — K] 


(36) 
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where K denotes the functional corresponding to the operator K. After parametrizing 
= [3oe^u^ and integrating Eq. (36) with respect to vr, we obtain the path-integral formula 
for the Massieu-Planck functional as 


= In / 


(37) 


with 


r ~ r 1 

dr / d^-^Xy/^N 

rf^o 


idr<p — 


i2N^ 

dr / d‘^~^x^/^C{(l),dp(l)]gf,y), 




(38) 


where N = NvPe^ = —nufd^/fdo, We dehne the thermal metric gpp and its 


inverse g^’^ as 


-N^ + NiN^ Ni 


9piv ~ 


r. 




-iV-2 N-^N^ 

Ni 773 / \ N-^N^ 7 ^^" - N-^N^N^ 

Here, Aj = jijN^ = This metric again has the form of the ADM metric: 

= —{Ndiy + ji-j^N'^di + dx''){N~^di + dx^), 


(39) 


(40) 


with dt = —idr. In Fig. 2, we show a schematic hgure of a locally thermalized state by com¬ 
paring it with that of the globally thermalized one. While the (uniform) thermal field theory 
is formulated under the flat spacetime as shown in Fig. 2(a), the locally thermalized held 
theory can be formulated under a curved spacetime background. The metric is determined 
by the thermodynamic parameters such as the temperature, and the huid four-velocity as in 
Eq. (39), and thus the imaginary-time radius manifestly depends on the spacetime as shown 
in Fig. 2(b). The line element, ds^ is not real because dt is imaginary, so that the action 
S[(j), A] is in general complex, which causes the sign problem in lattice simulations. This 
expression of the thermal metric does not explicitly depend on the choice of the original 
shift vector A*. 

The thermal metric is invariant under the imaginary-time translation, since the param¬ 
eters A“ do not depend on the imaginary time t, t = —ir. Furthermore, we also have local 
symmetry by the redehnition of the imaginary time. In order to demonstrate this symmetry, 
we rewrite the thermal metric from the ADM form to the Kaluza-Klein one as 


ds^ = —e‘^'^{dt + a-idx^Y + 'Jjjdx'^dx^, 


(41) 
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(a) Global thermal equilibrium 


(b) Loeal thermal equilibrium 



FIG. 2. Comparison between the global thermal equilibrium (a) and local thermal equilibrium 
states (b). 

where aj = 'j'jj = qq- + ujuj, and we used = —e^'^. In this 

parametrization, the square root of determinant of metric becomes \/—g = N^ = e^’-y/y. 
This parametrization of the Massieu-Planck functional was discussed in Ref. [28] . Following 
Ref. [28], we can easily see that this metric is invariant under the local transformation (the 
Kaluza-Klein gauge transformation), 

f ^ + x(s), 

(42) 

a-iix) a-iix) - dixix), 

where x(*) is an arbitrary function of the spatial coordinates. We note that qqj nonlinearly 
transforms under this transformation since yT does not change, so that 7 is not gauge 
invariant. This symmetry enables us to restrict possible terms that appear in the Massieu- 
Planck functional [28]. For example, aj appears in the Massieu-Planck functional only 
through the gauge invariant combination such as the field strength, j-q = — djai. 

In addition to the above symmetry associated with the imaginary time translation, the 
Massieu-Planck functional has the {d — l)-dimensional spatial diffeomorphism, x —)■ x'[x). 
This spatial diffeomorphism invariance also restricts possible terms that could appear in the 
Massieu-Planck functional. For example, 7 ' appears only in combination with i.e., 

= dTjiNe~^. In Sec. IV, we will write down the possible form of the Massieu- 
Planck functional within the derivative expansion using these symmetric properties. 

Although we only consider the neutral scalar field, the extension to a system with finite 
chemical potential is straightforward: We may replace the partial derivative d^- with the 
covariant one, = [dr — in which the additional term = z///9o is Kaluza-Klein 
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gauge invariant. Therefore, the symmetric properties of the thermal metric, which are 
discussed in this subsection, also hold for systems with hnite chemical potential. 


III. TIME EVOLUTION 

In the previous section, we considered the local thermodynamics on the hypersurface. 
Here, we discuss the time evolution of the expectation values of local operators. In a quantum 
held theory, the expectation value of a local operator is given by 

{d{x)) =tipod{x), (43) 

where po is the density operator at initial time. In particular, we consider the time evolution 
of hydrodynamic variables Ca{x). If the constitutive relation is obtained, i.e., if {Ja) is 
expressed as a functional of Ca or A“, its time-evolution equation (hydrodynamic equation) 
is given by the continuity equation = 0. To obtain the constitutive relation, it is 

useful to decompose {Ja) into nondissipative and dissipative parts, {Ja) = {Ja)\'^ + {^Ja) ■ 
The nondissipative part {Ja)];^ is obviously a functional of J'{x) and does not contain the 
information of the past state. On the other hand, we need the information of the past 
to evaluate {SJ^). The purpose of this section is to derive the self-consistent equation to 
determine {5J^). 

At a very early stage of time evolution, the system will be far from equilibrium in a state 
that cannot be characterized by only thermodynamic or hydrodynamic variables. In this 
stage, microscopic degrees of freedom play an important role to determine the time evolution 
of the system. In contrast, at later times, we expect the system to be characterized by the 
thermodynamic variables whose time evolution is governed by the hydrodynamic equations. 
In this paper, we assume that at the time the distribution function is given by a local 
Gibbs one, po = Plg[^o; although, in general, this is not exact but only approximate. As 
we will see below, once we assume this initial condition, the time-evolution equation can be 
rewritten as a compact form. 

In order to evaluate the expectation value of 6J^{x) at the point x^ G Sf for t > to, we 
decompose the density operator into the local Gibbs distribution on Sf and the other; 

p(fo) = exp(-^[to; A]) = exp(-^[t; A]-h S[t, to; A]), (44) 
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where S[t,to;A] = *5[t; A] — ^[to;A]. S[t, to;A] can be expressed by the divergence of the 
entropy cnrrent operator as 

[ dsds 

to 

In the last line, we used Eq. (30). The explicit form of is given in Eq. (32). 

We will treat S[t, to] A] as the perturbation term in the derivative expansion because 
is proportional to the derivatives of the parameters, V^A“. In order to expand p{io) with 
respect to S[t, to;A], we decompose the density operator as 


S[f,fo; A] 




ds / dE.NVJ^. 


(45) 


'to 


P{to) = PLG{t)U{t,'to), 


where U{t, to) is dehned as 


(46) 


U{t, to) = T,e^o dr±Ai,to-,\]^ (47) 

with fo;A] = e'^^^^’^^t,[i,io] Here, denotes r ordering. The expectation 

value of an operator 0{x) on E* is given by 


{6(x)) = {u6{x))A, (48) 

where {0{x))]:'^ = tr plg[^; A](9(x). If one takes O = U~^, Eq. (48) gives an identity 
corresponding to an integral fluctuation theorem, {U~^) = 1 [30]. 

Since Eq. (48) is the identity, it holds for any parameters A“. We need a condition to fix 
A“. Here we impose {ca{x)) = {ca{x))]:^ [31]; they are explicitly 

n^x{x){f^),{x)) =nf,{x){f^^{x))^^, (49) 

n^x{x){J^{x)) = n^(a:)(>(a:))p. (50) 

The parameters are determined by the entropy functional through Eq. (23). Equations (49) 
and (50) mean that the dissipative parts {SJ'^) are orthogonal to n^{x), i.e., n^{SJ'^) = 
-{Sca) = 0 . 

In order to consider the time evolution, we use the spatial projection operator introduced 
in Sec. HA, 


Pl^ = di; + with = -1, 


0 , Pj^n^ = 0 . 


(51) 
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Then, the derivative is written as 


Ti 

v„ = (-v’-n, + F;)V„ = -^Vr+ (52) 

where V* = Nv^V^ and V_l^ = By using this projection operator, S[t, A] reads 

as 


S[t, to] A] 





(V,A“)hc, + 


(53) 


We would like to eliminate the time derivative of parameters VsA“ from E[f, fo;A], which 
can be performed by using the continuity equation, V = 0- 
Since 5[t;A] does not depend on x\ i.e., V_L^^[t;A] = 0, A] = —{n^/N)diS[t] X], we 

can write the divergence of as 




LG 

t 



d^N{x') [ 

Jo 

dS,N{x')(V,\'‘(x'))(6c^{x),Sjt%x'))[, 


( 54 ) 


where (A, B)i is the local Gibbs version of the Kubo-Mori-Bogoliubov inner product, 

iA,B)i= f drie^^Ae-^^BAj^, (55) 

Jo 

which satishes linearity {aA + bB, C)i = a(A, G)f + b{B, C)i, Hermite symmetry (A, i?)| = 
{B, A)i, and positivity (A, A)* > 0; (A, A)i = 0 ^ A = 0. We used ((ca(a:))f SJ^ {x'))i = 
{ca{x))\^{5J}^{x'))\^ = 0 to obtain the last line in Eq. (54). Using Eq. (52), we hnd that 
= 0 leads to 


j dm6ca{x),5c,{x'))tVtX\x') 

+ J d^i6ca{x),5jax%N{P)V^yX\P) + N{x)V,{5XAA) = 0. 


(56) 


Multiplying Eq. (56) by the inverse of {Sca{x),6cb{x'))i, and integrating it with respect to 
the coordinates on the hypersurface, we obtain 


VfA“(x) = -j d^J d^{6c^{x),6cb{x^))A{Scb{x^),6j:{xA)iN{xAV^yX^{xA 

- [ dJ:'^{6ca{x),Scb{x'))AN{x')VMA^'))- 


(57) 
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Let us eliminate in T,[t, to; A]. For this purpose, it is convenient to introduce a projection 
operator V onto 6ca, 


VO = / dSt / dT!i5ca{x){5ca{x),5cb{x'))O[3ci,{x'),0)t. 


(58) 


This is the relativistic version of the projection operator used in Refs. [38, 39]. At thermal 
equilibrium, it reduces to the Mori projection operator [40]. We have 
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(S'- ( \ S' ( d\\x') 

{dCa[x)OCb[x))^ = 


SCa(x) ■ 

Using Eqs. (59) and (60) and the chain rule, we can rewrite Eq. (58) as 


Now, by using "P, we can eliminate VfA“ from S[t, t^o; A], and we obtain 


S[t,to; X] = -J dsJ dS,-iV[(Va;.A“)(l - V)5J,^ - 

ni /» 


LG 


(59) 

(60) 

(61) 


(62) 


^ - J da J dE,N 


- 5A“V„{5J-;) 


Here we introduced 50 = (1 — V)50, which enables us to remove the hydrodynamic modes 
from SO. In the second line, we replaced {SJ'^) by {SJ^) because the expectation value of 
the projected operator vanishes, {VO) = 0. We also dehned 


5A“(x) = f 


(63) 


5cb{x')' 

For later use, we perform the tensor decomposition for ST^'^. Since n^5T^^ = 0 and 
ny5T^’^ = 0, we can decompose 5T^^ as 5T^^ = h^^^Sp + where 

1 


6p = 






= P^P^STP^ 


hl^'^ 


(64) 

(65) 


Here we introduced and that satisfy h^^^hpy = Pp. 

As a result, S[t, Pq; A] reads 

±[t, to; X] = -j_'dsj dE,N [{h^’^VpP.)~6p + + (V±m^)5> - SX^Vp{Sjp) 


( 66 ) 
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where 


h 

V(^/3.) ^ ^(Vp/3. + V./3p) - (67) 

We note that V^{5j^) does not contain the explicit time derivative of the parameters 
because = (-A'“'ivV,-+Vi„)((5Ji‘> = (iV“‘(Vtn„) +Vi„)(i j"), where we used 

jr) = -(Vi7r„){i j:). 

Since (5A(x)) = (5A(a:)), onr goal is now to solve 

{5J^{x)) = Ae7o(68) 

St-[ 7,to;-^] contains {Sjj^^x)) as in Eq. (66), so that Eq. (68) becomes a self-consistent 
eqnation. As we discuss in the next section, {^jj^ix)) can be evaluated order by order in 
the derivative expansion with respect to the parameters. 

IV. DERIVATIVE EXPANSION AND HYDRODYNAMIC EQUATIONS 

In this section we perform the derivative expansion to derive relativistic hydrodynamic 
equations order by order. We also discuss the frame choice, which originates from an ambi¬ 
guity in the dehnition of the fluid four-velocity. 

A. Derivative expansion 

The expectation value of Ja{x) consists of the nondissipative and dissipative parts, 
{ja{x)) = + {SJ'^{x)). As will be shown in the following, the nondissipative 

part (j^aix))]:^ is obtained by differentiating the Massieu-Planck functional T with respect 
to i. The Massieu-Planck functional can be expanded as 

OO 

>I>[A]=5^>I>(">[A], (69) 

n =0 

where n denotes the order of spatial derivative 0(V([) As was discussed in Sec. II C, T[A] 
and therefore enjoy thermal Kaluza-Klein symmetry and spatial diffeomorphism 

^ On curved space, curvatures may appear in higher-derivative terms. For example, we identify the spatial 
curvature as the second-order derivative, because it is given by a commutator of the spatial covariant 
derivatives. 
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invariance. For parity symmetric theories, vanishes because we cannot construct a 

scalar with one spatial derivative such that it is invariant under the above symmetries. On 
the other hand, the higher-order terms are not forbidden by parity symmetry. The second¬ 
er higher-order hydrodynamics can contain nondissipative terms coming from them. The 
general expansion of the nondissipative parts based on these symmetries was discussed in 
Ref. [28]. 

The dissipative part {Sj'^{x)) can be expanded as 

OO 

(WW) = (70) 

m,n=0 

where the term labeled by (m, n) contains m temporal derivatives, Vf, and n spatial deriva¬ 
tives, V_L. In order to evaluate {6j'^{x))(^n,m), we expand the dissipative part {6j'^{x)) 
as 

{SJ^{x)) = (T^e^o 

= + f dT{TX(t,h)Sj!:{x))\° (71) 

+ \ j dr [ dT'{Trt^(i,h)tr:(t,io)SJ"{x))\° + --- . 

^ do do 

Here {5Ja{x))\^ vanishes by dehnition. Since to) contains the derivative of the param¬ 
eters, V_lA“, J^rii,to) is identified as of order V_l. We note that J2r(t,io) does not contain 
the temporal derivative of the parameters, VfA. This fact implies that the derivative ex¬ 
pansion starts from (^77^(a:))(o,i); i-e., for / > 0 vanishes. If one considers 

the nth order of {SJ'^{x)), one may expand Eq. (71) up to the nth order of S^(f,to)- All 
correlation functions with lower orders of J^rit,io) contribute to the nth order of {Sj'^{x)). 
For example, in addition to the third term in the second line of Eq. (71), the second term 
contributes to {5Ja{x))(o,2) through the derivative expansion of the correlation function 
{Tr^r(t,to)5JI^{x))\'^. In the following, we restrict ourselves to the zeroth and first-order 
hydrodynamic equations with parity symmetry. 

1. Zeroth order: Perfect fluid 

Let us consider the leading order of (fTfl) in the derivative expansion. We show that the 
energy-momentum tensor and the current have the form of a perfect fluid. In Sec. IIC, we 
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discussed that the Massieu-Planck functional is obtained from the path integral in curved 
spacetime, whose metric is invariant under the thermal Kaluza-Klein transformation. 
Thanks to the Kaluza-Klein gauge symmetry, does not contain aj. Furthermore, 

the spatial diffeomorphism invariance restricts the 7 dependence of to the form 

proportional to while it does not restrict the a dependence of \k^°)[A]. Then, we 

factorize as [28] 


r-/3o 


^(o)[A] = / dr 


= / dT^i 


(72) 


where jd' = (3 = /doe'^, and is the pressure of the perfect fluid as explicitly 

shown later. To obtain the second line, we used the relation /S-y/Y = /5'Y7 the fact that 
the parameters are independent of the imaginary time. 

Next, we consider the variation of ip with respect to t, which changes the hypersurface 
and We obtain 


dip = d{/3'p) = p^dfd^ + n'du — jd^pdn^ 
By using this relation, the time derivative of A] reads 

= f dEjN 


(73) 


1 

(v,/?Yp + /SYv,^)Yp + ^vyYp) 


= j d^iN 




(74) 


+ - p/?YV.nY) 


= / dSW 




Comparing Eq. (29) with Eq. (74), we obtain the expectation values of the energy-momentum 
tensor and the particle current by the local Gibbs distribution as 




{T^vi.x))t = S'pp +—{pu + n^p) 
= {e + p)u^u,y 6^p, 


= Y— = 7W^, 


(75) 

(76) 
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where e = = puU'"(3', and n = —n'(3^Uf^/(3'. Here we used that the 

energy-momentum tensor is symmetric under the change of the indices to derive the second 
line in Eq. (75). Equations (75) and (76) are nothing but the constitutive relations of the 
energy-momentum tensor and the particle current in a perfect fluid. 


2. First order: Navier-Stokes equations 

Let us consider the next leading order in the derivative expansion. We need not con¬ 
sider the derivative corrections coming from since vanishes for the parity- 

symmetric system. The hrst-order correction to the dissipative part comes from 

f drirtAt, t(i, (77) 

Jo 

We used the inner product Eq. (55) and the Hermite symmetry of the inner product. Then, 
the first-order corrections read 


(hf'^"(x))(o,i) ^ (78) 

(^>(a:))(o,i) ~ CSJ^{x),f:{t,io))i. (79) 

where ~ denotes an equality at the hrst order in derivatives. The right-hand side of 
Eqs. (78) and (79) also contain the higher-order contributions. In the hrst order in the 
derivative expansion, we can neglect (5A“V^(577^) in S(f, to) because {SJ'^) = 0(V) and 
thus V = O(V^). We can replace K in these inner products with P^(3^[x). We 
remark here that the dissipative corrections are orthogonal to by construction, and thus 
we do not need to employ or v'^ for the tensor decomposition. Therefore, we may decom¬ 
pose these inner products in Eqs. (78) and (79) by only using Two-point correlation 
functions with odd numbers of indices, such as {Sp{x), SJ^{x'))i, vanish. Furthermore, corre¬ 
lation functions with a single 5ti^^{x) also vanish since is traceless. In consequence. 
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we have 


{6p{x),t{t,to))i = 


{6n^%x),±it,io))t = 


(h>(a;),S(t,to))f = 


-j_ dt' j dT.pN\5p{x)Mx'))th^"{x')^M^') 

-j‘dt' j <iEpV{M'“'(x),«''-’(x')),-V(p/3,)(x') 

2p 


(3{x) 


h'^"(a;)h""(x)V(A)(x), 


dt' / dJ:pN'{U^^{x),Sr{x'))iV±^iy{x') 


'to 

K 


l3(x 


-V>{x), 


( 80 ) 


( 81 ) 


( 82 ) 


where we used 9^A“(x') ~ d^\°‘{x). Here the transport coefficients, (, p, and k, are the bulk 
viscosity, the shear viscosity, and the diffusion constant, respectively. They are given by the 
Kubo formulas: 


( =/3{x) / dt' / dJ:pN'{Sp{x'),Sp{x))i, 

J —oo J 

m 


V = 


K = 


{d+l){d-2)j_ 

P{x) 


dt' / d'Lt'N'{5%:^^''{x'),5%P'"{x))th^p{x)K„{x), 


d-l 


dt' / dJ:pN'{6J^{x'),6r{x))thp,{x), 


(83) 

(84) 

(85) 


where we replaced to by —oo, which can be justihed in the first order in the derivative 
expansion. We can now construct the constitutive relations up to the first order, which 
are given as Eqs. (5), and (6). Once we calculate the transport coefficients, (, p, k, and 
the pressure p{(3, u) from the microscopic theory, we have closed equations composed of 
the continuity equations. These are nothing but relativistic versions of the Navier-Stokes 
equations. We emphasize here that we derive them without choosing a frame such as the 
Landau-Lifshitz or Eckart frame. 


B. Choice of frame 

In relativistic hydrodynamics, we face the frame ambiguity, which stems from a way 
to define the fluid four-velocity. One useful frame is the Landau-Lifshitz frame, in which 
the energy flux of a fluid element vanishes at the rest frame of the fluid. Another is the 
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Eckart frame, in which the particle flux is absent. In our approach, the choice of and 
corresponds to the choice of frames. In this subsection, we show that by explicitly choosing 
and n^, our constitutive relations reproduce those in the Landau-Lifshitz and Eckart 
frames within the derivative expansion. 

1. Landau-Lifshitz frame 

The fluid four-velocity in the Landau-Lifshitz frame is dehned by the condition that in 
the local rest frame, the energy flux of a fluid element vanishes. Then, the energy and charge 
densities coincide with the local thermodynamic values. In other words, the Landau-Lifshitz 
frame is dehned by [1] 

{5f^'^{x))uLv{x) = 0 , {5J^{x))uL^,{x) = 0 , ( 86 ) 

where the subscript L denotes the Landau-Lifshitz frame. We can easily see that Eq. (86) 
is satished if we choose = u^. In this case, we have a familiar projection 

The constitutive relations up to hrst order in the derivative expansion 

read 


(f^"(x)) = (e + p)u^Lu'i + (87) 

(>(a;)) = n< - (88) 

where 

+ V^ULa) - 9 = V^<. (89) 

In this frame, we can explicitly write down the projected operators in Eqs. (83)-(85) as 

(91) 

~ /V /V 77 

SJ>^ = - h^'^5p,. (92) 

e + p 

To derive these equations, we used 

= (I) «+(!),« +0(Vx). (93) 
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n 


V5J^= / rfSf / d^5pp{x){5pp{x),5p,{x'))^\5p,{x'),5Jn-t = h^^''5p,^^ + 0{y^), 
J J e + p 


( 94 ) 


where 5e = —u^5p^, and we used the following relations [27]: 


J dSi(5pp(a;), 6pu{x'))i = ^hp^{e + p) + 0(V±), (95) 

J dE-tidp^ix), 5P{x'))t = + 0(Vx). (96) 


2. Eckart frame 


Next, we consider the Eckart frame. The fluid four-velocity for the Eckart frame is 
dehned by the condition that it is proportional to the particle current, i.e., u'P{x) = 
J^{x)/ J^(a;)J^(a;), where the subscript E denotes the Eckart frame, and J^{x) = 
(J^(x)) [9]. It is also required that the energy density is expressed as e = Up{T^'^{x))u^{x). 
In the hrst order in the derivative expansion, we may choose and as 


= 




K 






K 




Using = u^ + {K/{Pn))d^u -|- O(V^), we obtain 


(f ^'^(x)) = (e + p)u%ul + pg^^'^ + - Cdh^\ (98) 

{J^{x)) = nu'P, (99) 


where we dropped the terms of order V^. and 6 are obtained by replacing in Eq. (89) 
with u'p. The thermal conductivity which is absent in the Landau-Lifshitz frame reads 


q 


/4 _ 


e + p 
n(3 




( 100 ) 


We note that the shear and bulk viscous terms are the same as those of the Landau-Lifshitz 
frame. 

Although we do not have the charge diffusion in this frame, the expression of heat current 
is slightly different from the original Eckart one q'P, which is given by [9] 


g^ = -A(V((T + TViU^), 


( 101 ) 
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where A denotes the thermal conductivity of the fluid. The apparent difference is coming 
from whether we use the time derivative of the fluid four-velocity in order to construct 
the constitutive relations. Although we utilize the Mori projection operator to eliminate 
the time derivative of the parameters from the entropy production, we can reconstruct the 
constitutive relations by using the time derivative terms with the help of the equation of 
motion. In the first order, we can use the equation of motion for the perfect fluid, 

1 tlT 

( 102 ) 

T e + p 

in order to eliminate from Eq. (100). Then, we derive the constitutive relations in the 
original Eckart frame with A = ((e + pY/'n?)K. 

Obviously, in our formalism, the constitutive relations in the Landau-Lifshitz and Eckart 
frames are equivalent within the hrst order in the derivative expansion. These are related 
to each other by the redehnition of the fluid four-velocity, -H- -f (^/(/dn))in 
Eqs. (87) and (88). More generally, if we choose a frame such that + 0(V) and 

= M^-|-0(V), the constitutive relations in this frame are equivalent to those in the Landau 
frame within the hrst order in the derivative expansion. Namely, if is a functional of A“, 
the constitutive relations are unique and become those in the Landau-Lifshitz frame. We 
note that such a uniqueness was also discussed in Ref. [41] based on the Boltzmann equation. 


V. DISCUSSION 

In this paper, we have derived hydrodynamic equations from quantum held theory by 
assuming that the density operator has the form of the local Gibbs distribution at initial time. 
In particular, we have derived the hrst-order equations, that is, the relativistic version of 
the Navier-Stokes equation without a choice of frames such as the Landau-Lifshitz or Eckart 
frame. Our frame-independent analysis becomes important if the vorticity is nonzero; in this 
case, we cannot choose because the vorticity, constructed from and pU^, 

vanishes by the Frobenius theorem [32]. 

The real-time evolution in our formulation is schematically shown in Fig. 3. The density 
operator of the system at initial time to is assumed to have the form of the local Gibbs 
distribution. Then we expand the density operator at a later time t around the new local 
Gibbs distribution with the thermodynamic parameters A“(a;) at that time. In each time. 
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the local Gibbs distribution (the Massieu-Planck functional) can be expressed by using 
the imaginary-time path integral under the curved spacetime background Sf x S^, whose 
metric is given in Eq. (39). After a sufficiently long time, the system reaches the global 
thermal equilibrium with the uniform imaginary-time radius /3o- The local Gibbs distribution 
enables us to treat a nonequilibrium state beyond the real-time formalism [42] , in which the 
distribution is necessarily in the global equilibrium. However, in an early stage far from 
equilibrium, the density operator cannot be approximated by the local Gibbs distribution, 
and thus our formulation is no longer applicable. 



FIG. 3. Schematic figure of the real-time evolution in our formulation toward the global thermal 
equilibrium. 


As mentioned in the Introduction, our method is closely related to that presented by 
Sasa [30]. In fact, if we take = (—1,0) in the flat spacetime, they are equivalent. The 
difference is that our formalism is based on the Heisenberg picture, while that in Ref. [30] is 
based on the Schrodinger one; these are related to each other by the unitary transformation, 

;.LG,Sasa _ 

Pt — ^ Pt ^ 

There are several directions on future research based on this method: One is the gener¬ 
alization to a system with a quantum anomaly such as chiral fermions in which the matter 
couples to external gauge helds. This generalization is straightforward: We may replace the 
energy-momentum tensors and the particle current to those in background gauge helds. In 
this case, K[t] A] formally has the same form as before. The difference is that the currents 
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are no longer conserved, 


(103) 

V^> = (104) 

where F^^ is the held strength of the external gauge held, the antisymmetric tensor, 
and Cano the anomaly coefficient. Using Eqs. (103) and (104) instead of Eqs. (1) and (2), 
the divergence of the entropy operator reads as = —(Vu/3^)ST^^ — 6J^{x)fp, where 

fp = VpV + P'^Fpi,. The term Fp^^Fp^ coming from the anomaly cancels out in 

the divergence of the entropy current operator. Therefore, the anomaly does not directly 
contribute to the dissipative part of the currents, which is consistent with the observation 
in the entropy-production method [43] and the generating-functional method [28, 29]. 

Another direction is an application to second-order hydrodynamic equations. There are 
several works derivingthese equations from microscopic theories [12-19], which are based on 
the Boltzmann equation. In contrast to them, our method is applicable to strongly coupled 
systems. In general, all possible terms respecting symmetries appear in the derivative ex¬ 
pansion, whose coefficients depend on details of the system. Our method gives the Kubo 
formulas for these coefficients. We may obtain Kubo formulas different from those in the 
analyses based on the Boltzmann equation. We leave these interesting applications for future 
work. 
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Appendix A: Derivation of Eq. (30) 


Let us, here derive Eq. (30). Noting that the volume element can be written as = 
dfixy/^5(t — i{x))d^t{x) = —d^x^/^^fJ_6{t — i{x)), we write 


j = - j d'^xy/^d^e{t-t{x))f>^{x) 

= j d‘^x^/^e{t-t{x))\/i,f^{x), 


(Al) 


where we used the integral by part, and assumed that f^{x) vanishes at the boundary. The 
derivative of Eq. (Al) with respect to t leads to Eq. (30), 


di j dE,J^{x)=dj j d‘^x^e{i-i{x))VJ^{x) 
= j dJ:tN{x)V,nx), 


(A2) 


where we used dSf = d^Xy/—g6{t — t{x))N 
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